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Abstract 

The characteristics of the partial nuclear muon capture with mas- 
sive left-handed Dirac neutrino and relativistic component of the muon 
wave function have been derived. The multipole amplitudes are given 
as a function of neutrino mass parameter and reduced nuclear matrix 
elements which are modified by the small component of the muon wave 
function. As an example, the capture rate, asymmetry and polarization 
of recoil nuclei are investigated in terms of these multipole amplitudes. 
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1 Introduction 



Neutrino mass problem is a topic of many theoretical and experimental works. In recent 
years a considerable effort has been devoted to find massive neutrino in the different 
processes. It has been extensively looked for in neutrino oscillations, beta decay, double 
beta decay (3(3 (especially neutrinoless beta decay (3(3(0v)), charged lepton conversion as 
well as in cosmology and astrophysics (see e.g. [I- 4]). 

Massive neutrinos are naturally predicted beyond the standard model of electroweak 
interaction, e.g. by " see-savJ 1 mechanism. One can imagine also the gravitationally in- 
duced neutrino masses [5]. 

Experimental data of (3(3 and (3(3{Qv) decay give limits on the mass of Majorana neutrino 
and right-handed currents [6]. The solar neutrino puzzle can be explained in frame of 
MS W mechanism [7,8], or by long- wavelength vacuum oscillations [9]. 
As for muon capture process, the possibility to verify the existence of massive muon neu- 
trino was proposed in ref. [10]. The calculations of nuclear muon capture observables 
with massive neutrino were undertaken in refs. [11,12]. So far, however, the calculations 
were performed without the small relativistic corrections which a priori could be com- 
peted with neutrino mass effect. We present here the results of our calculations of the 
observables of nuclear muon capture, taking into account both the neutrino mass effect 
and relativistic component of the muon wave function. 



2 Effective Hamiltonian and multipole expansion 

As earlier, we use the method of multipole expansion [13, 14] for the effective Hamiltonian 
in the nuclear muon capture on the lepton-hadron level, i.e. Fermi theory of four-fermion 
point interaction. The calculation is performed with relativistic muon and massive left- 
handed Dirac neutrino. 

We take the effective /^-capture Hamiltonian in the standard current x current form 

1 1 

<(l,^,z\H f ,\-m,z>=GJ x lx, (1) 

where G = Gm ^ s 6c is the coupling constant, J\ is the hadron current operator, q is 
vector of neutrino momentum, and 

h = *V7a(1 + 75 )^ (2) 

is matrix element of the modified lepton current with the massive left-handed Dirac 
neutrino. With z, unit vector along z axis, we denote the direction for lepton spin 
quantization. 

The Dirac bispinor is commonly used for the neutrino wave function, 
1 . 1 ... /-Ci(<T • q) 



^Uf, =< r | q, -v, z >= —j= exp{+iq • r} y ^ M; ) <g> Xv , (3) 
where c and C\ are energy-mass coefficients, 



H 1 -* (4) 
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We consider the nuclear partial transition in the muon capture process; 

»- + (A,Z)j i ^(A,Z-l) J/ + v li - (5) 
With the method of [13], we get the multipole expansion; 

<-q,J / M / ;qi/ i ,q|//J J,M i; >= gr^LT?" ^ ^J^V^)" ( 6 ) 

The multipole amplitudes (h, fi = ±|, | Jj — J/ |< L < Ji + Jf ) are given by linear 
combinations of the coupling constants and the reduced matrix elements in the nuclear, 
or in general, hadron space. With q, we denote that lepton spin is quantized now along 
the direction of neutrino momentum. 

We find, that neutrino energy-mass coefficients factorize from the multipole ampli- 
tudes as follows; 



1 L — 



c — 2hc 1 



nhfl 



(7) 



It should be noted that still are matrix elements of lepton current, where the diagonal 
terms Tj^ h are amplitudes with no spin-flip of the lepton in the weak current, while T^~ h 



are the amplitudes where the lepton spin is flipped in course of interaction. These 
amplitudes essentially differ for the neutrino helicities h — +\ and h — —\- We have, 
like in [13], the following expressions for the multipole amplitudes in terms of the 
reduced matrix elements in hadron space; 



rphh 

1 L 



(8) 



[OLL,J 4 ] hh + i[lLL,J] hh +4 



2L + 1 



[lL-lL,3] hh +4 



I L + l 
2L + 1 



[!L + lL,J] hh , 



T, 



h-h 



(9) 



[OLL, J 4 ] h -h ~ i[lLL, J] h - h +\ 



I L + l 
2L + 1 



[lL-lL,J] h _ h +4 



2L+ 1 



[lL + lL,J] h _„. 



Here, [OLL, J^\h±h an d [lL'L, 3]h±h are notations for the reduced nuclear matrix elements, 
in the sense of the Wigner-Eckart theorem. 

For the relativistic muon in the Is atomic state (K — shell), we take the Dirac 
bispinor for the particle bound in the radial electrostatic potential; 



-iF(r)xim(r) 
G(r)x-im(r) 



(10) 



with the customary definition of spherical spinors [15], 



X fcm (r) = Y l (r) ® xi/s 



Jm 



(J = I - -sign/c). 



(11) 



G(r) and F(r) are named as large and small radial components, respectively. 
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Next, these matrix elements are given explicitly in the product form of two terms, 
where the former is bilinear combination of the radial functions which emerge from the 
matrix elements of the leptonic current due to the neutrino plane wave and Is muon 
bound state, and the latter is tensor operator that in the standard way couples the 
neutrino spherical harmonics and hadron current operators: 

[0LL,J 4 } hh = (12) 
^ < / II G{G(r)j L (qr) - 2hF(r)[^— [ j L ^(qr) - ^jjL + i(qr)]}J,Y L (r) \\ i >, 

[0LL,J 4 ] fc _ fc = (13) 



J2 JL(L + 1) 

JJ<f\\ G(-2h)F(r) y 2L + 1 [yL-i(gr) + j L+1 (qr)} J 4 Y L (r) \\ i >, 

[lLL,3] hh = (14) 

V2 JUL + 1) . . 

< f II GF (r V 2L + 1 [h-Mr) +J L+1 (qr)}[J®Y L (r)] L \\ i >, 

[lLL,J] h _ h = (15) 
^ < / || G{G(r)j L (qr) - F(r)[A^-i(?r) - L + i(?r)]}[J ® Y L (i)] L \\ i >, 

[lL-lL,3] hh = (16) 

^ < / || G(-2h)[G(r)j L ^(qr) + 2hF(r)j L (qr)][J ® Y L -\r)] L \\ i >, 
J f 

[LL-lL,3] h _ h = (17) 

^ < f || G(-2h)[G(r)j L - 1 (qr) + F(r)j L (qr)][J®Y L - 1 (r)} L \\ i >, 

[lL + lL,3} hh = (18) 
^ < f || G(-2/*)[G(r)j L+1 (gr) - 2hF(r)j L (qr)][J ® F L+1 (r)] L || i >, 

[lL + lL,J] ft _ fc = (19) 
^ < / || G(-2h)[-G(r)j L+1 (qr) + F{r) 3L {qr)}[3 ®Y L+ \v)] L \\ i >. 

jv(qf) are the neutrino spherical Bessel functions. Tensor couplings of the hadron current 
spatial components J and spherical harmonics from the neutrino plane wave are as 
follows, 

iJ®Y L \i)] L M = ^Wyf(r). (20) 

jitn \ ' / 
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3 Relativistic muon and massive neutrino 



Question now is for the multipole amplitudes to compare the leading relativistic muon 
(RLMU) effects with bare massive neutrino (NEMA) contribution. To answer this, 

at least qualitatively in 0** > l Wf transitions, we have included RLMU effects into 
the formulae for multipole amplitudes calculated within the framework of the impulse 
approximation for the hadron current. Then, we generalize Fujii-Primakoff form factors 
to; 



G k y = gy(l-2h&), 

G\ = gA + 2hgv{l + n P -n n )^, (21) 



G\-G% = -2hg A + (g A -gp + 25i h gp%. 



2M 



It is convenient for our present purposes to use point nucleus solution (Coulomb poten- 
tial). Then the well known relation holds; 



F(r) = _^i_2 G (r)~-^G(r), (22) 



where 7 = yl — (aZ) 2 . a is a fine structure constant and Z - atomic number. We have 
considered contributions due to F(r) together with derivative G'(r). 
Reduced matrix elements are written within the short notation; 

{lL'LJttr} = ^< f\\j^G{r i )j l {qr i )[a i ®Y L '{v i )] L r-\\i > (23) 
Jf i=i 

and additionally {1L 1 L, jyc} = [1L'L,<t], if the labels in Bessel function and spherical 
harmonics are the same. 

Below, the multipole amplitudes in 7Vi — > l w f transitions with the leading order RLMU 
contributions, ~ are given; 

= (G A -G h P )^([m,a}+V2[121,a})+g A [0U,a-i^} (24) 
+ ^-9A^l({101,j l( T} - ^{121,^0"}) , 

iT^ h (^f) = -2hG h A ([101,0-] -a/2[121,<t]) - g v [Ul,ij^] (25) 
- 2h%-g A ({101,^0-} - ^{121,^}) , 



Thh I \ 
1 \^i=-^f) 



G$[0U, 1] + 2hg v] jl ([101, i±] + v/2[121, i-^]) (26) 
+ %- ■ - (2^ y {011, j l} - 9aV2{111, j Q <r}) , 
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Tt%^f) = -G h Ani,a}+2hg v ^[101,i^}-V2[121,i^}j (27) 

+ — f2^ y ^-{011,jol} + ^2{lll,joo-}j • 

The largest RLMU contributions are given here. These are only from the small 
radial component of the muon function F(r). All reduced matrix elements with Bessel 
functions of rank L — 2 and L — 3 have been omitted as they can not contribute for 
multipolarity L — 1. 

Estimating NEMA contributions we perform with Lommel integral for spherical 
Bessel functions for (q < E) in case of massive , and (q = E) massless neutrinos; 



oo 



7T 



r(s±i) i 



h(q) = Uqr)dr = • • " (28) 



o 



2 r(^) g 



With this, the difference between j v (qr) and j v (Er) in radial integral with nuclear wave 
functions can be evaluated; 

IM = -^rUE) * (l + l ^§] I V (E) (29) 
c Ci \ 2 q 1 J 

Multipole amplitudes change approximately by a global factor (1 + |^f). We take 
the upper limit m u c 2 < 160 keV from [16], and q ~ lOOMeV. The magnitude of 
= 1.28 x 10~ 6 measures the change of the amplitudes. Coulomb coupling parameter 
aZ for e.g. 12 C gives \aZ ~ 2.16 x 10 -2 . This suggests RLMU effect to the multipole 
amplitudes by four orders of magnitude stronger than NEMA contribution. 



4 Observables in muon capture 

In order to obtain the capture rate and correlation characteristics in muon capture, we 
use the polarization density matrix for the final state of the system, 

< -q, J f M f ; q, l -v\p s \ -q, J f M' f ; q, V >= (30) 
= 2(2j] + 1) ^ < ~ q ' 3 ^ M ^\ V I H » I 3 ^\ m >< \ m I l + 2P^-s M | l -m > 

Mi,m,m 

< -q, J f M' f ; q, V | Hp \ JM; l -m >* , 

where direction of z is the spin quantization axis for both nucleus and leptons. is the 
residual muon polarization on the K-orbit and s M = \a denotes the muon spin operator. 
Averaging over the initial states and summing over the final ones, which are not measured 
in the experiment, we obtain the formulae for nuclear muon capture rate, A c ~ Trpj = 

JdPftj^f dqTrpf, recoil polarization vector P/ = Tr (J/P/) /JfTrpj, and the recoil 
asymmetry W(0) = Trpf/Tr~pf. 
The formula for muon capture rate is; 
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here coCiNfi is a phase space volume, with conventional factor for m u = 0; 

-l 



1 + 



yjm* + (Am P Y 



(32) 



q is absolute value of the neutrino momentum, and CqC\ = [1 + (m^/g) 2 ]^ 1 / 2 . A is the 
nucleus mass number, while and m p are muon and proton masses, respectively. All 
the multipole amplitudes come into the coefficient, 



a, = e (i n h i 2 + 1 n- h i 2 ) 



ih-h l 2 ^ 



(33) 



Then, we define the weights u>h', 

( 1 — 2hcoCi\ 
= 



A, 



(34) 



These weights correspond respectively to the contributions of the neutrinos with h — +| 
and h — — \ to the observables. In the lowest order approximation for the energy-mass 
coefficients, eq. (4), (with m v /q <C 1) we get uo- ~ 1, and oo + ~ 0. 

Here we give the recoil polarization formula in the Jj = — > J/ = 1 transitions; 



ft 



(35) 



where 



P^/i) = [-2/iA ft + (A - tfeS/J (P M • q)] q + tfe/^P^ + 2WmB h (P (J x q). (36) 

The structure functions A ft and are expressed in terms of the multipole amplitudes 
ratio; 



hh 



h-h 



-x h exp [i(j> h ], 



(37) 



2 + xV 



B h 



2x h exp [i(p h ] 
2 + x 2 h 



For the longitudinal polarization, P L = j- J rfq < Pj • q >, we get 

P L = J2^h2hA h . 



(38) 



(39) 



Next we consider angular asymmetry of recoil with respect to the direction of the muon 
polarization; 



W(0) = l + 5> h 2/ia h q.P,,, 

h 



(40) 



7 



where P M cos 9 = — q ■ P M . We express the asymmetry coefficients ah as the functions of 
x h and find their relation to recoil polarization coefficient A h ; 



a h = §— ^ = - 1. (41) 

A -|- Xi 



With the results of [12] for the neutrino polarization, we present simple example of the 
link between CP-violation and massiveness of the Dirac neutrino. We have calculated 
in accordance with the definition < s„ >= Tr (s u pf) /Trpf\ 

< s„ >= [a + (c- Reb){P^ • q)] q + RebP^ + Imbq x P M , (42) 

where structure functions are given by, 

a = Y,huj h , c = -]-^2uj h a h , (43) 



2 ' h 



m u 1 , 

^i — ^i T? 



1 + c Ci \ / 1 - C Ci \ 

) X + ) X " 



(44) 



(Since then, we denote h = sign h and ji = si^n /x). 

So that, observing in muon capture < • q x >~ Imb ^ would indicate both 
CP-violation and massive neutrino. For the longitudinal neutrino polarization we have; 

< Sj, ■ q >= a — cP^ cos 9, (45) 

which gives for m u = the recoil asymmetry; < Sj, ■ q >= — |Wo(0). Averaging this over 
P M , we get the helicity for the left-handed massless neutrino; < • q > av = —\- 



5 Discussion and conclusions 

In the model with massive left-handed neutrino, [12], the existence of the four indepen- 
dent multipole amplitudes was obtained (doubled in comparison to the massless case) in 
the nonrelativistic approximation of the muon wave function. As follows from our con- 
sideration the calculation with accurate muon wave function does not change the number 
of multipole amplitudes, which is still four. So the form of the kinematic characteristics 
of the muon capture expressed by multipole amplitudes remains the same. However, 
the numerical predictions will change, as the number of nuclear matrix elements in the 
multipole amplitudes is larger owing to relativistic muon contribution. There appear 
additional nuclear matrix elements [ILL, J] hh and [OLL, J^h-h- Moreover, the structure 
of other nuclear matrix elements is also different due to additional terms proportional to 
F(r) and G'(r). 

Considering nuclear observables, we obtain strict results in terms of the multipole 
amplitudes. In addition to standard term, there is another one with the amplitudes for 
positive helicity of neutrino. This term multiplies by the weight uj + , which is proportional 

2 

to In these observables the multipole amplitudes with opposite neutrino helicities 
do not interfere. 

Additionally , NEMA effects contribute to neutrino Bessel functions in the multipole 
amplitudes and very weakly to the phase space volume. From these, the effect to Bessel 
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function shifts the multipole amplitudes by a term proportional to and doubles its 
full value in the capture rate. 

On the contrary, the ratio of the multipole amplitudes, 

T hh/ T h-h^ ( 46 ) 

weakly depends on NEMA and RLMU effects, which practically disappear in allowed 
transitions, i.e. for 7Tj = 7rj. In conclusion, such observables as recoil polarization, recoil 
asymmetry etc., to a high degree are independent of these effects. 

Qualitative analysis suggests relativistic muon effects larger by the four orders of 
magnitude than massive neutrino effects in light nuclei. In progress, there is our analysis 
on specific nuclei for quantitative results. Situation is different in muon capture reactions 
with three-body final states, like e.g.; 

(A, Z) 3i + iT - (A - 1, Z - l) Jf + n + v fl , (47) 

There are experiments performed and planned, 

H~ + 2 H -> z/^ + 2n, 

+ 3 He -> i/p + 2 H + n, 

to measure high energy-momentum transfer to hadrons, which automatically gives low 
energy neutrino, m u /E v — > 1. Eventual analysis of the data with the massive neutrino 
theory could be possible. 
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